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The parametric representation is given to the multisoliton solution of the Camassa- 
Holm equation. It has a simple structure expressed in terms of determinants. The proof 
of the solution is carried out by an elementary theory of determinants. The large time 
asymptotic of the solution is derived with the formula for the phase shift. The latter 
reveals a new feature when compared with the one for the typical soliton solutions. The 
peakon limit of the phase shift is also considered, showing that it reproduces the known 
result. 
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§1. Introduction 

In this paper, we report some new results associated with the multisoliton solution 
of the Camassa-Holm (CH) equation 1 ) 

u t + 2k 2 u x - u txx + Suu x = 2u x u xx + uu xxx . (1) 

Here, u = u(x, t), k is a positive parameter and the subscripts t and x appended to u denote 
partial differentiation. Originally, this equation has been found in a purely mathematical 
search of recursion operators connected with the integrable partial differential equations. 2 -* 
Recently, equation (1) has attracted considerable interest since it has been derived as 
a model equation for shallow- water waves. 1 ) In addition, the equation has been shown 
to be completely integrable. With this as a turning point, a large number of works 
have been devoted to studying the mathematical structure of the equation. A recent 
paper describes a short history and relevant references concerning the CH and related 
equation. 3 ) Almost all the works have been focuced on the case k = for which the CH 
equation exhibits peakon solutions which are represented by piecewise analytic functions 
and whose dynamics are now well understood. 4 ) However, when k ^ 0, several new features 
appear in solutions. In particular, solutions recover analytic nature, but they are expressed 
by a parametric form like u = u(y,t),x = x(y,t) where y is a new coordinate variable. 
A difficult technical problem is to find the inverse mapping x = x(y,t). To date, this 
problem is resolved only for particular cases. 3 ' 5 ' 6 ) In this respect, we remark that an 
approach based on the inverse scattering transform method (1ST) provides an explicit 
form of the inverse mapping in terms of Wronskian determinants. 7 ) Nevertheless, the 
general iV-soliton formula is not available yet. 

The main purpose in this paper is to present a complete description of the general 
AT-soliton solution in the form of parametric representation. We show that the solution 
exhibits a simple structure expressed in terms of two determinants. In §2, we shortly 
summarize the derivation of a system of equation (associated CH equation) equivalent to 
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equation (1). In §3, we describe the main result in which a parametric representation is 
given to the iV-soliton solution of the CH equation. The proof of the solution is carried 
out by an elementary theory of determinants. In §4, we explore the asymptotic form of the 
iV-soliton solution for large time and show that the interaction process of solitons reveals 
a new feature when compared with the typical one for the Korteweg-de Vries equation. In 
§5, we summarize the results and discuss some open problems related to the CH equation. 

§2. Associated CH equation 

It has been demonstrated that the CH equation can be recast into a more tractable 
form by means of an appropriate coordinate transformation. 6 ) Here, we give a short sum- 
mary of the results. Introducing the new variable r in accordance with the relation 

r 2 = u - u xx + k 2 , (2) 

the CH equation (1) can be put into the form 

r t + {ur) x = 0, (3) 

where the boundary condition for r is r(±oo, t) = k. Then, we define the coordinate 
transformation (x,t) — > (y, t') by 

dy = rdx — urdt, dt' = dt. (4) 

In the following analysis, we use the time variable t in place of t' by virtue of the second 
relation in (4). Transforming (3) by means of (4), it becomes 

r t + r 2 u y = 0, (5) 

and u is expressed in terms of r as 

u = r 2 — r(ln r) ty — k 2 . (6) 
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We term the system of equtions (5) and (6) the associated CH equation. 6 ) If we substitute 
(6) into (5), we obtain an alternative but more convenient form in the following analysis, 
which is 

Qt = r y , (7) 

where 

q = yt-\(t) + \{h~h)- (8) 

By eliminating the variable r from (7) and (8), we can see that Q evolves according to 
the following nonlinear wave equation 3 ) 

Qt + 2k 3 Q v + 4 K 2 QQ t + 2K 2 Q y d~ 1 Q t - K 2 Q tyy = 0, (9) 

where d~ l = — f°° dy is an integral operator. An important observation is that equation 
(9) can be identified with a model equation for shallow-water waves. 3 ) This fact enables 
us to obtain the iV-soliton solution of (9) in the (y, t) coordinate system. To complete the 
solution, however, one must revert to the original (x, t) coordinate system via the inverse 
mapping 

x y = 1 x t = u(y, t). (10) 
r(y,t) 

The most difficult ingredient in the analysis is to integrate (10) for the iV-soliton solution. 

§3. Parametric representation of the iV-soliton solution 

Now, the main result in this paper is summarized as follows: The iV-soliton solution 
of the CH equation (1) can be written in a parmetric representation 

u(y,f)= ( ln f) f > (11) 

x = ^ + ln fy + d. (12) 
« h 

Here, fi = fi(y,t) and f2 = f2(y,t) have the determinantal expressions 

/ 1 = |G|,/ 2 = |tf|, (13) 
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with the N x N marices G and H whose elements are given respectively by 

G = {9ij),9ij = -^-e%- + = 1,2,.., N), (14) 

Pi -q% Pi- qj 

H = (htj), hij = -^-e%- + = 1, 2, .., N), (15) 

Pi-Qi Pi- <l.j 



where 



& = (Pi - Qi)y + £ f- - - ) * + &o, (i = 1, 2, .., AT), (16) 



2 VP; ?■ 



Pi = ^(1 + Kfei), 9i = ^(1 - (i = 1, 2, .., iV), (17) 

Sij is Kronecker's delta and the parameter d in (12) is an integration constant. By virtue 
of the parametrization (17), the iV-soliton solution is characterized completely by the 2N 
parameters hi and ^q, (i = 1, 2, N) . In terms of the parameters fcj, the phase variable 
£i of the ith soliton may be written in the form 

2k 3 

1 ~K 2 k 



i 



& = ki[y- q ^Tf < - Vio ) > (* = 1. 2, iV), (18) 



where we have put £jo = — ^ij/io- 

Let us now outline the proof of (11) and (12) in which two bilinear identities (23) and 
(38) below will play an essential role. First, we write the iV-soliton solution of equation 
(9) in a determinantal form 3 ' 8 ' 9 ) 

g = -2(ln/) w ,/=|F|, (19) 

where F is an N x N matrix with elements 

e& 1 

F = (fa), fij = + — — , = 1,2, .., AT). (20) 

Pi Qi Pi Qj 

Substituting (19) into (7) and integrating the resultant equation by y under the boundary 
condition r — > k, \y \ — > oo, we obtain 

r = « - 2(ln /) tv . (21) 
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A crucial observation in the present analysis is that r can be represented in terms of /, fi 
and fi as 

(22) 



r = k- 



f 



2 ' 



It follows from (21) and (22) that 



Kfif2 = Kf 2 -2(ff ty -f t f y ). 



(23) 



This is a bilinear identity among determinants. We note that analogous relations have been 
studied in the direct proof of the multiperiodic solutions of the Benjamin-Ono and nonlocal 
nonlinear Schrodinger equations while employing an elementary theory of determinants. 10 - ) 
For later convenience, we first introduce some notations as well as formulas for determi- 
nants and then describe the main result. Matrices and cofactors associated with any 
N x N matrix A = (ay) are defined as follows 

I an 

A(cii;bi) = 



ai 



A(^cii, hi] Ci, di) — 



/ an 



ajvi 
ai 
V h 



CIin 



A , 



a N / 
a 1N ci di ^ 

ttNN Cn dN 


b N / 

d 2 \A\ 



(24) 



ij,kl 



(25) 



(26) 



da ik daji 

Here, A t j and Ai^ ki are the first and second cofactors, respectively. The following formulas 
are used frequently in the present analysis: 

an ... a 1N xi 



(inn xn 
Vn z 



N 

= \ A\z — ^ ] AijXiyj, 



|A(a,,6,; Cl ,^)||A| = \A( ai ; c i )||A(6 < ; d t )\ - \A( ai ; dJWAfc; c;)|, 
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(27) 



(28) 



N N 

J^ifi + gMjAij =J2(fi + 9i)\A\, (29) 

i,j=l i=l 
N 

E (/»• + 9s)a rs A is A rj = (fi + gj)Aij\A\. (30) 

r,s=l 

Formula (28) is Jacobi's identity and formulas (29) and (30) follow from the expansion 
formulas for determinants, Ylk=i a ikAjk = 8ij\A\,J2k=i a kiAkj = Sij\A\. 

Let us now prove (23). Using the rule for the differentiation of determinant and 
formula (29) with A = F, fi=Pi and = —q^ we obtain 

N N 

\F\ y = -J2 F ij + J2(Pi-Qi)\F\. (31) 

i,j=l i=l 

Similarly, we find from (29) with A = F, fi = p~ x and gi = —q~ X that 



N T? N / 1 1 \ 

r K J ' K I ) \F\ (32) 



\t 



If we use formula (30) with A = F, f r = p r and g s = —q s as well as Jacobi's identity of 
the form \F\Fujj = FuFij — FijFu, we can derive the formula for the y derivative of Fij 

{ N N I f N r r 

Fa* = - E ^ + E(p' - ^ §i + E n^r " {Pi " Qj)Fij - (33) 

[ l,s=l 2 = 1 J l ' s=1 

Differentiating (32) by y and inserting (31) and (32), we obtain the expression for \F\ ty . 
With (31) and (32), this result is substituted in the right-hand side of (23) to obtain the 
relation 

N N N _ 

~ 2(ffty ~ ftfy) = ^ - « £ ^ E — + E (34) 

Further simplication is possible by applying (27) to (34). This gives 

nf - 2{ff ty - f t f y ) = k(\F\ + iF^p-'DQFl - \Fiq~ 1 ; (35) 
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On the other hand, using (27) together with the basic formulas for determinants, we can 
show that 

N 

fi = \G\ = n(Pi/fc)(l*1 + m^Pi^l), (36) 
i=i 

N 

f2 = \H\ = H(q j /p j )(\F\-\F( q - 1 ;l)\). (37) 
i=i 

The identity (23) follows immediately from (35), (36) and (37). 
The identity below also plays an important role: 

r 2 ./,./••> /.•(./,./••>., ./,.,/,.). (38) 

which we shall now show. First, it follows from (36) and (37) that 

f - hh = WW 1 ; 1)1 - impi 1 )]) + W^W^- i)|. (39) 

Formulas similar to (31) now take the form 

N N 

\G\ y = ]J( Pj /q j )\F( qi -,p- 1 )\ + J2(Pi-Qi)\G\, (40) 

j=l i=l 
N N 

\H\ y = H(q j /p j )\F( q - 1 ;p i )\ + J2(Pi-Qi)\H\- (41) 

j=l i=l 

Then, we calculate the quantity P = f 2 — fif 2 + n(fi, y f2 — fifi, y )- Substitution of (36), 
(37), (40) and (41) into P yields 

P = \F\(\F(qT' 1 ; 1)| - K\F{q~ 1 ;p i )\ - \F{\;p^)\ + K\F{q i ;p~ 1 )\) 

Let Pi be the sum of terms involving |F| and P 2 be the rest. Owing to the basic formula 
for deteminant, a\F{a,i m ,bj)\+(3\F{ai m ,Ci)\ = \F{ai\ abi+(3ci) \ and the relation 1 — npi = nqi 
which follows directly from (17), P\ reduces to 

P L = K \F\(\F(qr 1 ;qJ\-\F( p r 1 m )\). (43) 
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Using Jacobi's identity and (17), we can show that 

P 2 = «|F||F(p i ,l;pr 1 ,«r 1 )l- 



(44) 



After a few manipulations, we finally arrive at the relation 



P = P 1 + P 2 = 



k\F 



(\F{q % - P% )\-\F{ Pl - qi )\), 



(45) 



i=l piQi 



where F = (fij) is an N x N matrix with elements 



fij =Pi(Qifij + 1), (i,3 = 1,2, N). 



(46) 



Thanks to (17) and (20), we see that F is a symmetric matrix and hence \F(qi;pi)\ 
|P(Pii9i)l) implying that P = 0. Thus, we complete the proof of (38). 



The relations (22) and (38) immediately lead to our main result. Indeed, it follows 
from (22) and (38) that 



and substitute (47) into (48) . Integrating the resultant equation in y under the boundary 
condition u — > 0, \y\ — > oo, we obtain the expression (11) for u(y,t). In view of (47), a 
system of equations (10) which determine the inverse mapping can be integrated. This 
gives rise to (12). Note that the parameter d is independent of t as confirmed by (48) and 
the second equation in (10). The expressions (11) and (12) give a complete description of 
the parametric representation for the iV-soliton solution of the CH equation. In the process 
of constructing soliton solutions, we have noticed that it is better to use the determinants 
fi and f 2 in place of the determinant /. It should be pointed out, however, that fi 
and f 2 exhibit the same functional form as / except phase factors. Te see this, let / = 




(47) 



We rewrite equation (5) in the form 




(48) 
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...,£ N ). Then, it follows from (14), (15) and (20) that h = f{^ 1 + ln{p 1 /q 1 ), ...,6v + 
\n(p N /q N )) and f 2 = /(£i + ln(q 1 /p 1 ) : ...,£ N + \n(q N /p N )). 



§4. Asymptotic behavior of the iV-soliton solution 

A new feature of the iV-soliton solution appears in the interaction process of solitons. 
Here, we address on this problem. The procedure for investigating the asymptotic behavior 
of the solution can now be performed straightforwardly. The core part of the calculation 
is to evaluate fi and f 2 by utilizing the formula for the Cauchy determinant 

(Pi ~ Pj)(qi ~ Qj) 



d(m, m + 1, n) 



Pi ~ Qi 



n 



m<i<j<n 



(Pi - Qj)(Qi - Pj) 



(49) 



,Pi Qj J m <i,j<7 

To this end, we order the magnitude of the velocity of each soliton in the (x, t) coordinate 
system as c\ > c 2 > ... > c/v where 

2k 2 



Ci = 



(z = l,2,...,iV). 



(50) 



We take the limit t — > — oo with the phase variable & of the ith soliton being fixed. Since 
then other phase variables behave like ^1,^2, ~^ +00, £h-2, ...,£jv — > —00, fi 

has the leading-order asymptotic of the form 

\ / N 



fi ~ exp 



i-i 



X 



^e ?l rf(i + 1, % + 2, AT) + i + 1, JV)^ 
By invoking (49) and (17), we obtain 

d(i+ l,i + 2,...,N) 



(51) 



d(i,i + l,...,N) ^{ki-kj 



(52) 



Substitution of (52) into (51) gives 



fi ~ exp 



i=i 



AT 



n- n— Wi+i,...,iv)(^!- , + i), 



(p Ai 



Ql 



\S = 1 



(53) 
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where 



,(-) 



N 



7* ' = £ ln 

j=i+l 



k^ kj 



,(i = l,2,...,N). 



ki ~\~ kj 

Similarly, in the limit of t — > — oo, f 2 has the asymptotic form 



(54) 



f 2 ~ exp 



i-l 

i=i 



n^)(n p -^)<M +1 ,...,.)(|e- : -- +1 ). 



(55) 



It turns out from (11), (53) and (55) that u is represented by a superposition of N solitons 



i=i 

where Ui is a one-soliton solution given by 3 ' 5 ' 6 ) 



(56) 



*2iKiCi k^ 



1 + k 2 A;2 + (1 - k 2 A;2) cosh £ 
In the same limit, the mapping relation (12) becomes 



(57) 



6 (1 - /cfeiJe^-Ti ' +1 + /^ 

x — Cit — Xio ~ — — h ln 



' vA '< ' (1 + /cA^e**-^ + 1-^ 
1 + nki 



-ln 



+ d, (z = l,2,..,iV), 



(58) 



1 — Avfe 

where = Uio/K- In the limit of t — > +oo, on the other hand, the expressions corre- 
sponding to (56), (54) and (58) are given respectively by 

N 



«~£«i(&-7, (+) ) 



<+> = !> 



i=i 



1 A ,(i=l,2,...,N), 



&i "I - kj 



(59) 



(60) 



& ^ (l-Kkj)e^ +) +1 + Kkj 
x - Cit - x i0 ~ — + ln — } 

Kk i (1 + Kki)e^-^i +1- nki 



N 



E, / 1 + k/c 



K/kn 



In 



1 + ftfcj 
1 Kikj 



+ d,(i= 1,2,. .,N). 



(61) 
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Let Aj be the phase shift of the zth soliton in the (x,t) coordinate system. This quantity 
can be evaluated simply with an appropriate use of (58) and (61). The result is 

A-^-l-Vg.n^) 2 - E ■n(l±g) , ,( i -l,2,..,«). W 

1 j=l v 3 7 j=i+l V 3 7 

The phase shift consists of two contributions. The first term on the right-hand side of (62) 
comes from u{y, t) and the rest terms from x(y, t). Note that the first term is the same as 
the phase shift arising from the interaction of N solitons for the KdV equation. However, 
due to the mapping (12), additional terms appear as indicated by (62). The case N = 2 
in (62) recovers a formula already derived by the 1ST. 11 ) 

It is also an interesting problem to investigate the characteristics of the iV-soliton 
solution in the limit of k — > 0. Here, we are concerned only with the limiting form of 
the phase shift. We find that the appropriate limiting procedure can be performed if one 
puts Kki = 1 — a and takes the limit ti — > while keeping Ci(i = 1,2, ...,N). At the 
same time, the conditions = ejCj(i,j = 1,2, ...,N) must be imposed by virtue of the 
relation between Cj and ki (see (50)). The formula (62) then reduces to 

i-l , \ 2 N , x 2 

A ' = -E ln M^ + E ln U^" .(' = 1.2. ••.")• «53) 

j=l \ 1 3 ' j=i+l ^ J7 

In the special case of N = 2, (63) coincides with a formula presented in ref. 4). 
§5. Discussion 

We have obtained a simple parametric representation for the iV-soliton solution of 
the CH equation. The solution is expressed compactly by the two determinants fi and 
ji. This finding is of essential importance. If one calculates the solution u from (16) 
and (21) in terms of a single variable /, then resulting expression becomes a complicated 
fashion and it would not be of practical use in investigating the structure of the solution. 
Introduction of the variables fi and fi also leads to an analytical form (12) for the inverse 
mapping. The explicit form of the iV-soliton solution makes it possible to construct other 
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class of solutions. For instance, the rational soliton solutions may be obtained from it by 
taking appropriate long wave limits hi — > 0(i = 1, 2, N). 11 ^ 

Recently, the CH equation has been generalized to a two-dimensional version by ap- 
plying an asymptotic expansion method to a system of water-wave equations. 12 ) It is an 
interesting problem to investigating its integrability. The method developed in this pa- 
per may be useful in constructing multisoliton solutions, if they exist. Furthermore, the 
Degasperis-Procesi (DP) equation is a current interest in soliton theory. 13 ' 14 ) Although 
the DP equation has a form similar to the CH equation, its mathematical structure is 
quite different from that of the CH equation. 14 > Quite recently, we have succeeded in ob- 
taining the multisoliton solution of the DP equation by means of a reduction procedure for 
the multisoliton solution of the Kadomtsev-Petviashvili equation. 15 ) The solution can be 
written in a parametric form analogous to the corresponding solution of the CH equation. 
However, the simple expressions like (11) and (12) are not at hand yet for the general 
iV-soliton solution. This problem is currently being investigated. 
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